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I .  INTRODUCTION 

The  objective  of  this  paper  is  to  evaluate  the  relative  importance 
of  shock  heating  and  plastic  work  in  the  initiation  of  explosives  by 
projectile  impact.  There  is  some  confusion  about  this  issue  due  to  t!te 
large  variety  of  impact  tests  which  are  conducted  on  explosives.  To 
eliminate  some  of  this  confusion  and  to  delineate  the  situation  which  we 
are  considering,  we  will  first  describe  some  of  the  tests  which  have 
been  used  and  what  is  known  about  the  initiation  mechanism  in  each. 

A number  of  different  phenomena  fall  under  the  broad  topic,  "impact 
initiation  of  explosives".  The  term  impact  initiation  is  frequently 
used  to  describe  drop  weight  tests.  These  tests  are  distinguished  by 
having  a thin  sample  of  explosive  sandwiched  between  a large  drop  weight 
(projectile)  and  a rigid  anvil.  The  velocities  and  pressures  required 
for  initiation  in  such  experiments  are  very  low  (generally  less  than 
5 meters/sec  and  10  kb  respectively).  In  this  type  of  test,  Heavens  and 
1 2 

Fields  and  Afanasevand  Bobolev  have  demonstrated  that  initiation 
occurs  as  the  result  of  adiabatic  deformation  of  the  thin  samples 
(Figure  1) . Their  work  indicates  that  plastic  and/or  viscous  flow  play 
an  essential  role  in  the  initiation  process.  Mechanical  shock  is 
unimportant  under  these  conditions. 

In  a second  type  of  experiment  a flying  plate  impacts  an  explosive 
charge  which  is  sandwiched,  at  the  time  of  impact,  between  the  flying 
plate  and  a rigid  anvil.  In  this  ty])e  of  experiment  the  diameter  of  the 
flying  plate  is  large  compared  to  the  thickness  of  the  charge,  but  the 
charge  dimensions  are  much  greater  than  in  the  drop  weight  tests.  Such 

3 

experiments  have  been  performed  by  Napadensky  and  analyzed  numerically 

4 

by  Kot,  et  al  . In  the  experiments,  two  results  were  observed.  Above  a 
certain  velocity  detonation  occurs  within  a few  microseconds  as  a result 
of  the  incident  or  reflected  shock.  Below  this  velocity  explosions. 


1.  S.  N.  Heavens  and  J.  E.  Fields,  Proceedings  of  the  Royal  Society, 

A338,  77,  1974, 

2.  G.  T.  Afwiasev  and  V.  K.  Bobolev,  "Initiation  of  Solid  Explosives 
by  Impact,"  U.S.  Department  of  Cormerce,  Springfield,  VA,  1971. 

3.  H.  S.  Napadensky,  "Initiation  of  Explosives  by  Low  Velocity  Impact, " 
Fourth  Symposium  on  Detonation,  U.S.  Naval  Ordnance  Laboratory, 

White  Oak,  MD,  1965. 

4.  C.  A.  Kot,  H.  A.  Napadensky,  Y.  A.  Shikari,  and  A.  H.  Wiedemann, 

"A  Numerical  Study  of  Impact  Phenomena  in  Explosives  and  Propellants, " 
Proceedings  of  the  International  Conference  on  Computational  Methods 
in  Nonlinear  Mechanics,  University  of  Texas,  Austin,  TX,  1974. 
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} Figure  1,  Drop  Weight  Test  of  a Layer  of  Granular 

^ PETN  Viewed  in  Transmitted  Light. 

Interframe  Time  is  5.5  fxs. 

J (Photograph  Taken  From  Ref.  1 j 


but  not  detonations  may  occur  after  the  explosive  has  been  extruded 
laterally.  The  shock  initiation  process  has  been  studied  by  many  inves- 
investig- r.ors  and  is  qualitatively  understood.  Ignition  occurs  at  hot 
spots  which  are  formed  in  the  explosive  behind  the  shock  and  the  resultant 
energy  release  accelerates  the  shock  towards  detonation.  From  a micro- 
scopic point  of  view,  plastic  and  viscous  flows  in  the  neighborhood  of 
voids  are  necessary  preconditions  for  shock  initiation,  but  from  a 
macroscopic  point  of  view  deformation  is  not  important.  The  initiation 
of  low  order  explosions  is  not  so  well  understood,  but  it  seems  probable 
that  macroscopic  deformation  is  involved  in  a crucial  way.  Kot,  in  his 
numerical  analysis  of  the  experiments,  observed  that  material  strength 
had  an  insignificant  effect  on  the  temperatures  generated  by  the  shock 
in  the  explosive,  lie  concluded  that  energy  dissipation  due  to  distortion 
was  of  little  importance.  However,  his  calculations  did  not  extend 
to  times  large  enough  for  significant  macroscopic  distortion  to  occur, 
and  his  conclusion  does  not  apply  to  initiation  of  low  order  explosions. 

A third  experimental  geometry  involves  impacts  by  projectiles  which 
are  small  in  all  dimensions  compared  to  the  dimensions  of  the  receptor 
charge.  Bullet  impact  tests  are  typically  of  this  type.  Such  tests  may 
be  performed  with  either  confined  or  unconfined  explosives.  Broim  and 
5 6 

Whitbreau  and  Dewey  and  Slade  performed  and  analyzed  this  type  of  test 
for  charges  which  were  unconfined  or  confined  only  on  the  impacted 
surface.  With  such  charges  one  observes  either  that  the  charge  detonates 
or  is  shattered  and  the  result  is  recorded  as  no  reaction.  Dewey  and 
Slade,  and  Brown  and  Whitbread,  demonstrated  that  in  this  situation  the 
initiation  occurs  as  a result  of  the  shock  generated  at  impact  and 
macroscopic  material  deformations  are  not  important.  However,  when 
fully  confined  charges  are  impacted,  explosions  which  are  not  detonations 
can  occur,  and  these  reactions  can  occur  with  impact  velocities  which 

are  much  below  those  required  for  detonation^.  The  times  required  for 
the  occurrence  of  these  low  order  reactions  are  long  (several  hundred 
microseconds) , and  it  is  possible  that  material  heating  due  to  macro- 
scopic deformation  of  the  charge  may  play  some  role  in  the  initiation 
process. 

In  this  paper  we  focus  our  attention  on  the  latter  type  of 
experiment  where  the  charge  is  fully  confined  and  the  charge  dimensions 
are  large  compared  to  those  of  the  projectile.  We  are  particularly 
interested  in  the  mechanism  leading  to  low  order  explosions  as  opposed 
to  detonations.  At  the  onset  of  this  work,  we  expected  that  plastic 


5.  S,  Brown  and  E.  Whitbread,  Les  Ondes  de  Detonation,  2961. 

6.  J.  Dewey  and  D.  Slade,  BRL  Report  2021,  295?. 

7.  R.  Frey,  G.  Melani,  M.  Chawla,  and  J.  Trimble,  "Initiation  of 
Violent  Reactions  by  Frooectile  Impact, " Sixth  Symposium  on 
Detonation,  San  Diego,  CA,  2976. 


work  might  bs  an  important  mechanism  in  those  situations  where  the 
thickness  of  the  confinement  on  tho  impacted  surface  is  large  compared 
to  the  radius  of  the  projectile.  In  this  situation  the  shock  becomes 
weak  before  reaching  the  explosive.  We  were  interested  not  only  in  a 
comparison  of  shock  and  deformational  heating,  but  also  in  the  trends  of 
deformational  heating  with  projectile  radius  and  velocity.  In  this 
paper  we  have  considered  only  elastic-perfectly  plastic  materials. 
Viscous  effects  will  be  considered  later. 


II.  BRIEF  DESCRIPTION  OF  NUMERICAL  TECHNIQUE  FOR  EVALUATING 
THE  TEMPERATURE  RISE  IN  EXPLOSIVE 


A.  The  HELP  Code 


The  princi’'  .1  analytic  tool  employed  in  this  numerical  study  was  the 

g 

HELP  code,  which  is  a finite  difference,  multimaterial  Eulerian  computer 
program  dealing  with  compressible  fluids  and  elastic-perfectly  plastic 
flows  in  two  space  dimensions  and  time.  Although  the  code  is  basically 
Eulerian,  free  surfaces  and  material  interfaces  are  treated  in  a Lagrangian 
fashion  throughout  the  computational  grid.  No  material  diffusion  is 
permitted  across  these  discrete  interfaces.  In  the  past,  the  HELP  code 
has  been  used  to  solve  a variety  of  problems,  in  the  area  of  fragmenting 
munitions,  hypervelocity  impact,  and  shaped  charge  jet  formation. 

B.  Material  Model 


The  material  model  employed  in  HELP  consists  of  an  equation  of 
state,  a deviatoric  constitutive  relationship,  a yield  criterion  and  a 
failure  criterion  for  each  of  the  material  packages  employed  in  the 
problem.  The  materials  of  interest  here  are  steel  for  the  projectile 
and  the  casing  and  unreacted  solid  Comp  B3.  Table  I provides  values  of 
some  of  the  parameters  employed  in  the  calculation.  The  main  fV?atures 
of  HELP  arc  discussed  briefly  in  the  following  subsections.  The  CGS 
system  of  units  is  employed  in  HELP  code.  Therefore,  no  attempt  is  made 
here  to  convert  the  input/output  of  HELP  to  any  other  system  of  units. 

1.  Equation  of  State.  For  the  unreacted  explosive.  Comp  B3,  a 
Mie-GrUneisen  equation  of  state  was  employed.  This  equation  is  fully 
described  in  Appendix  A and  assumes  the  GrUneisen  parameter,  to  be  a 
constant.  The  form  of  the  equation  is 

P = (An  + B|i^  + Cvh  + PpE  (1) 


8.  L.  J.  Hageman,^  D.  E.  Wilkins^  H.  T.  Sedgwick  and  J.  L.  l/addell^ 
"HELP:  A Multi-Material  Eulerian  Program  for  Compressible  Fluid 
and  Elastic-Plastic  Flows  in  Two  Space  dimensions  and  Time^  " 
Systems,  Science  & Software  Report  No.  SSS-R-? 5-2654,  1975. 
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Table  I.  Elastic  Modulus,  Yield  Strength  and  Failure  Criterion 


where  P is  pressure,  E is  specific  internal  energy,  \i  = p/p^  - 1.  In 

order  to  ensure  that  the  pressure  is  tensile  for  expanded  states, 

(y  < 0),  equation  (1)  was  modified  to  read 

P = (Ay  - By^  + Cy^)  + PpE  (2) 

In  equations  (1)  and  (2) 

in  2 

A = 13.5  X 10  (dyne/cm  ) 

1 n 2 

B = 9.5  X 10  ^ (dyne/cm  ) 

in  2 

C = 100.6  X 10^^  (dyhe/cm  ) 

r = 0.947 


For  steel,  a Tillotson  equation  of  state  was  employed  in  order  to 
describe  the  condensed,  the  expanded  and  the  transitional  states.  For 
the  condensed  states,  i.e.  for  y > 0 or  for  any  cold  states,  E < E^, 
the  equation  has  the  form 


P = P = [a  + ] Ep  + Ay  + By^ 


+ 1 


(4) 


For  expanded  hot  states,  i.e.  for  y < 0 and  E > E',  the  equation  of 
state  has  the  form 


aEp  + 


bEp 


E n 
0 


+ 1 


+ Aye 


1) 

] 


(5) 


In  Equations  (4)  and  (^5)  n = p/p^  = y + 1,  A smooth  transition  between 

the  condensed  and  the  expanded  states  is  obtained  from  an  interpolated 
form  valid  for  E^  < E < E^  and  y < 0.  This  transition  region  equation 

of  state  has  the  following  form 


P = P^, 


(E  - E^)  Pg  . (E'  - E)  P^ 
^s-^s 


(6) 


In  equations  (1)  through  (6),  A,  B,  C,  a,  B,  a,  b,  E^,  E'  and  p^  are 

constants  for  the  particular  material.  The  values  of  these  constants 
for  steel  are  listed  in  Table  II. 
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Table  TI.  Equation  of  State  Constants  for  Steel 


a = .5 
b = 1.5 

12  / 2 
A = 1.28  X 10  dyne/cm 

B = 1.05  X 10^^  dyne/cm^ 

E = 0.095  X 10^'^  erg/g 
o 

a = 5 
B = 5 


00 

u 

o 

CL 

g/cm^ 

E = 2.44  X 10^® 
s 

erg/g 

E"  = 10.2  X 10^° 

erg/g 
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2.  Elastic-Plastic  Constitutive  Relation  and  Yield  Criterion.  The 
deviatoric  stress  increments,  ds..  are  determined  by  using  the  elasticity 
relationship 


2Gde . . 
11 


C7) 


where  G is  the  rigidity  modulus  and  de^^  are 
strain.  When  any  increment  of  stress  causes 


increments  of  deviatoric 
the  Von-Mises  yield  condition 


S..S..  < 2Y 
ij  11  - 


(8) 


to  be  violated,  each  stress  component  is  proportionately  reduced  to 
bring  the  stress  state  normally  back  to  the  yield  surface,  A variable 
yield  strength  incorporating  work  hardening  and  thermal  softening  is 
defined  to  be 


Y = (Yq  + YjM  + Y2P^)  (1  - |-)  (9) 

m 

In  equation  (9),  Y^  and  Y2  are  the  coefficients  of  work  hardening  terms 

and  E is  the  melt  energy, 
m 

3,  Failure  Criterion.  A material  is  said  to  have  failed  in  HELP, 
if  the  material  compression  falls  below  a critical  value  indicated  by 

9 

(p/p  )min*  The  critical  material  compression  is  defined  by 


(£_)  = . i + 1 

Po  MIN 


(10) 


where  S is  the  material  spall  threshold  and  k is  its  bulk  modulus.  All 
the  stresses  in  a cell  are  zeroed  out  if  its  compression  ratio,  p/p^  is 

less  than  (p/p  )„tm*  The  values  of  failure  criteria  for  the  two  materials 
0 MIN 

employed  in  the  code  are  listed  in  Table  I. 
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^ C.  Calculation  of  Plastic  Work 

I 

f In  order  to  discuss  the  calculation  of  plastic  work  in  any  material 

H package,  it  will  be  necessary  to  review  the  important  computational 

I steps  employed  in  the  HELP  code.  A complete  theoretical  description  of 

the  HELP  code  is  given  in  Reference  8.  For  each  computational  cell  of 


9.  R.  T.  Sedgwiak^  J.  L.  Wafiidellj  M.  Bakev  and  J.  M.  Walshs  "Studies 
of  Erosion  and  Impact  Resistance  of  ABM  Materials^  " BRL  Contract 
Report  No.  301,  1978. 
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the  fiulerian  grid  in  HELP,  the  conservation  equations  determine  mass, 
momentum  and  total  energy  at  any  time  cycle.  These  equations  are 


1^  /y  p dV  = -/g  p u^  n^  dS  (11) 

/y  p u.  dV  = /g  o. . n.  dS  - /g  p u.  u.  n.  dS  (12) 

Iy  /y  P E-p  dV  = /g  u^  n^  dS  - /g  D u^  E,p  n^  dS  (13) 

where 

Et  = 7 “i  Ui  + £ (14) 

is  the  specific  total  energy,  E is  the  specific  internal  energy, 

are  the  stresses,  u^  are  the  velocities  and  p is  the  density.  The  volume 

and  surface  integrals  in  above  equations  refer  to  integration  over  the 
volume  and  surface  area  of  the  cells  in  the  grid.  It  is  convenient  to 
express  the  total  stresses  as  made  up  of  two* parts;  hydrostatic  and 
dcviatoric,  i.e. 


(15) 


Employing  the  finite-differencing  scheme  along  with  equation  (15),  one 
obtains  from  equations  (11)  through  (13) 


Am  = - At  /„  p u.  n.  dS 

^ XI 


(16) 


A(m  Uj) 


At  /o  s..  n.  dS  - At  /o  P n.  dS 
S xj  1 S j 

At  /e  (p  u.  u.)  n.  dS 

O X J X 


A(m  E™)  = At  /q  s..  n.  u.  dS  - At  P u.  n.  dS 

1 ^3JXJ  o XX 


At  /g  (p  u^  E.J,)  n^  dS 


(17) 


(18) 


The  three  terms  on  the  right  respectively  represent  increments  on  the 
cell  surfaces  due  to  the  stress  deviatoric  forces,  the  hydrostatic 
pressure  forces  and  the  contribution  due  to  transport.  These  increments 
in  mass,  momenta  and  energy  are  calculated  in  three  distinct  phases  of 
the  computation.  These  phases  are 
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• S PHASE  where  the  effects  of  material 

strength  are  considered 

• H PHASE  where  the  effects  of  pressure 

are  taken  into  account 


• T PHASE  where  the  effects  of  material 
transports  are  studied 


In  the  absence  of  material  strength,  no  plastic  work  will  be  done  on  the 
material.  The  plastic  work  is  the  contribution  to  the  internal  energy 
from  the  S phase.  As  seen  from  equations  (16)  and  (17),  the  total 
energy  is  calculated  from  equation 


A (m  E„)  = At  / s. . u.  n.  dS 
s ''  r 13  3 1 

while  the  momentum  is  calculated  from 


(19) 


A 


s 


(m  u.)  = At  / s. . n.  dS 
r ij  1 


(20) 


Starting  from  equations  (19)  and  (20),  the  contribution  to  plastic  v/ork 
for  each  cycle  is  easily  obtained.  The  plastic  work  increment  of  any 
coll  i in  the  S phase  during  a time  cycle  n is  given  by 


A AC^  A 

A = A E.j,_  - A (j  u^  u^) 


(21) 


Because  of  the  material  flow,  the  plastic  work  and,  consequently,  the 
temperature  rise  could  be  grossly  in  error  if  care  is  not  taken  to 
follow  the  material  motion.  In  this  study,  a feature  of  the  plugging 
version  of  HELP  was  employed.  In  the  plugging  version  four  Lagrangian, 
massless  tracer  particles  are  placed  in  each  cell  of  a subgrid  of  the 
prviblem.  The  subgrid  is  chosen  to  contain  the  portion  of  the  material 
package  that  undergoes  most  deformation.  In  time,  the  tracers  move  with 
the  velocity  of  the  material  and  migrate  to  different  cells.  An  area- 
weighted  average  is  used  to  update  the  plastic  work  of  the  material 
associated  with  each  tracer.  Consider,  for  example  a point  0 in  Figure  2. 
The  plastic  work,  W^,  associated  with  point  0 at  time-cycle  n is  defined 
as 


Z A.  AW" 

1 -ill 

w"  = ^ 

0 4 


i=l 


A. 

1 


(22) 
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where  is  the  overlap  area  of  the  fictitious  cell  centered  at  the 
point  0 with  cell  i,  AW?  is  the  specific  plastic  work  increment  associated 
|i  with  center  of  cell  i during  cycle  n. 

j In  the  plugging  version,  once  the  plastic  work  associated  with  one 

i or  more  points  in  the  region  of  space  spanned  by  the  tracer  exceeds  a 

I predetermined  cutoff  value,  PLIVMIN,  the  plug  is  allowed  to  propagate. 

; Propagation  of  the  plug  was  inhibited  in  this  study  by  arbitrarily 

j setting  the  plastic  work  cutoff  value  too  high. 

! Since  the  plastic  deformation  of  the  material  is  irreversible,  it 

• is  safe  to  assume  that  all  of  the  plastic  work  is  converted  into  heat, 

j The  temperature  rise  due  to  plastic  work  is  then  easily  calculated,  for 

i a cell  centered  around  the  point  0,  from  the  equation 


i 


i 

I 

H 

i 


i 

t 

K 

% 


^.n  , 
0 ■' 


J C 


(23) 


In  equation  (23) , J is  the  mechanical  equivalent  of  heat  and  is  the 

specific  heat.  It  should  be  noted  that  each  cycle  makes  a positive 
contribution  to  equation  (23).  The  plastic  temperature  will  cease  to 
rise  whenever  the  material  fails  in  tension  or  melts.  In  both  cases 
shear  stresses  as  well  as  the  hydrostatic  pressures  are  zeroed  out.  The 
temperature  will  also  cease  to  rise  if  the  material  stops  deforming. 

D.  Calculation  of  Shock  Temperature 

yvhen  a strong  compression  wave  is  transmitted  into  metal,  the 
internal  energy  and  therefore  the  temperature  of  the  material  increases. 
It  is  assumed  that  material  strength  plays  a negligible  role  at  this 
stage  since  the  shock  pressures  are  usually  several  hundred  times  the 
yield  strength  of  the  material.  Moreover,  it  is  assumed  that  in  a time 

_7 

of  10  seconds  or  less,  thermodynamic  equilibrium  is  established  for 
the  material  behind  the  shock  front,  so  that  the  state  of  material 
behind  the  shock  front  can  be  described  by  an  equilibrium  temperature. 

The  Walsh-Christian  technique^®  used  here  to  calculate  the  shock  tempera- 
ture of  steel  and  explosive  is  described  in  detail  in  Appendix  C of  this 
report.  The  sliock  temperature  is  easily  calculated  from  the  expression 


T(V„)  = Tq  e 


b(Vo  - V„)  -bV„  V„  bV 
+ e / I ^ 

HUG 


dV 


(24) 


i 

}: 


10.  J.  M.  Walsh  and  R.  H.  Christian,  Rhys  Rev,  9^ 


2544,  1955. 
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where 


and 


b = r/v, 


0 


(25) 


f(V) 


dP 

dV 


(Vr 


(26) 


The  pressure  and  the  derivative  of  pressure  in  equation  (26)  are  readily 
evaluated  from  the  equation  of  state  of  the  material  under  consideration 
and  the  resulting  integral  in  equation  (24)  is  computed  numerically. 

After  attaining  the  peak  temperature,  the  material  expands  adiabatically, 
without  exchange  of  energy  with  the  outside  world.  Temperature  during 
this  process  of  expansion  therefore  drops  and  is  given  by 

b(V„  - V) 

T = Tj^  e " ' (27) 

where  (T^,  Vj,p  are  the  temperature  and  volume  of  the  shock  heated 

material,  while  (T,  V)  are  the  temperature  and  volume  of  the  adiabatically 
cooled  material.  In  order  to  calculate  the  residual  shock  temperatures 
of  the  materials,  equations  (24)  and  (27)  are  to  be  employed  with 
appropriate  code  values  for  the  volumes,  V^,  V^j  and  V. 


III.  RESULTS 

This  section  describes  the  results  of  a series  of  numerical  experi- 
ments that  were  performed  using  the  HELP  code.  The  quantities  of 
interest  are  the  interfacial  pressure,  shock  pressure  and  plastic  work 
temperature.  Since  the  explosive  in  this  study  was  always  confined  in  a 
steel  box,  it  was  therefore  possible  to  calculate  all  of  the  above 
quantities  for  the  explosive  as  well  as  the  steel. 

However,  it  must  be  stressed  that  the  numerical  values  obtained 
using  HELP  code  may  not  be  completely  dependable  until  a thorough 
analysis  of  computational  errors  introduced  by  finite  zone-size,  mis- 
matched-lmpedance,  artificial  viscosity,  etc.  is  made.  Also,  the  code 
cannot  predict  the  temperature  of  hot  spots  which  might  be  formed  by 
intergranular  friction  or  other  mechanisms.  The  purpose  of  this  study 
has  been  to  arrive  at  some  useful  conclusions  regarding  the  trends  as  a 
result  of  variation  of  impact  parameters  such  as  the  projectile  diameter, 
the  projectile  velocity  and  the  casing  thickness.  The  material  description 
employed  here  is  approximate,  at  best.  The  results,  therefore,  may  or  may 
not  match  the  experiments.  Each  calculation  reported  here  required  a core 
of  120  K decimal  -“'id  computing  time  in  excess  of  12  hours  on  a Univac  1108. 

A.  Typical  Case 

The  typical  set-up  of  the  problem  consisted  of  a cylinder  of 
unreacted  Comp  B3  explosive,  which  was  5 cm  in  radius  and  5 cm  high. 


The  explosive  wa<5  completely  enclosed  in  a 1/2  cm  thick  steel  casing. 

The  explosive/steel  assembly  was  impacted  by  a projectile  which  was  2 cm 
in  length  and  1 cm  in  diameter.  The  impact  velocity  was  1 km/sec.  The 
impact  took  place  along  the  axis  of  symmetry  of  the  problem,  thus  it  was 
sufficient  to  study  only  half  of  the  total  geometry.  The  computational 
grid  consisted  of  28  x 48  cells.  In  the  z-direction  the  cell  size  was 
0.2  cm  while  in  the  r-direction  the  cell  size  changed  progressively  from 
0.1  cm  to  0.5  cm.  The  cell  size  in  the  regime  of  greatest  deformation 
was  0.1  X 0.2  cm.  Figure  3 depicts  the  initial  projectile/explosive/casing 
configuration  emplcyed  in  this  calculation.  This  figure  also  shows  the 
four  tracers  per  cell  which  are  used  to  monitor  the  plastic  work  of  a 
cell.  The  tracer  particles  were  placed  only  in  a portion  of  the  entire 
grid.  Figure  4 shows  the  deformatio.  of  the  projectile  and  target  at 
60  lis  after  the  impact.  Figure  5 shows  the  compression  of  the  entire 
grid  at  t = 60  ps.  At  this  time  almost  all  of  the  cells  in  the  grid 
have  acquired  normal  or  less  than  normal  densities  indicating  extensive 
material  failure.  Figure  6,  which  is  a plot  of  pressure  in  the  projectile 
and  target,  indicates  that  most  of  the  high  pressure  regions  in  the  grid 
are  tensile.  Figures  7 and  8 show  the  radial  and  axial  velocities  at 
t = 60  ps.  At  this  time  the  proj<jctile  has  lost  most  of  its  velocity 
and  the  crater  lip  formed  in  the  target  is  moving  in  a direction  opposite 
to  that  of  the  projectile.  Figure  9 is  a pictorial  map  of  the  specific 
internal  energy.  Figures  10  and  11  are  the  maps  of  plastic  temperatures 
of  the  grid.  The  temperature  maximas  of  all  the  tracers  that  have 
migrated  into  a given  cell  are  plotted  in  Figure  10.  For  pure  cells, 

Figure  10  shows  either  the  maximum  temperature  of  the  steel  or  the 
explosive.  The  temperatures,  in  general,  are  lower  in  explosive  than  in 
steel,  mainly  because  the  specific  heat  of  the  explosive  is  roughly 
three  times  that  of  the  steel  and  because  the  explosive  yield  strength 
is  an  order  of  magnitude  lower  than  the  steel  yield  strength.  In  the 
mixed  cells,  containing  both  steel  and  explosive,  the  maximum  temperature 
will  correspond  to  a steel  tracer.  Figure  11  shows  only  the  temperature 
maxima  of  the  explosive;  the  steel  temperatures  are  completely  blanked 
out.  The  highest  temperatures  are  found  always  to  belong  to  the  tracers 
which  were  initially  placed  at  the  casing  interface.  In  60  ps,  most  of 
these  tracers  have  moved  away  from  their  initial  locations. 

Figure  12  is  a plot  of  the  highest  temperature  rise  versus  time  for 
both  the  steel  casing  and  the  undetonated  explosive.  For  each  material, 
a high  peak  in  temperature  rise  corresponding  to  its  shock  pressure  is 
obtained  in  a short  time  after  the  impact.  The  shock  temperature 
acquired  its  highest  value  in  the  steel  a little  earlier  than  in  the 
explosive  but  drops  quite  drastically  shortly  thereafter  in  both  materials. 
At  about  10  ps  each  of  the  materials  begins  to  deform  plastically  and 
starts  to  heat  up.  The  temperature  at  the  hottest  point  in  the  flow  is 
plotted  in  Figure  12.  The  highest  temperature  at  late  times  may  be  in  a 
portion  of  material  which  was  cooler  at  early  times.  In  other  words, 
the  temperatures  sketched  in  Figure  12  may  belong  to  different  portions 
of  material  at  different  times.  The  rise  in  plastic  temperatures  levels 
off  at  62°C  for  the  explosive  and  at  190°C  for  the  steel.  The  following 
observations  can  be  made  with  respect  to  Figure  12. 
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C Figure  3.  The  Projectile  / Explosive  / Casing  Configuration  at  the 

^ Beginning  of  the  Impact  for  the  Typical  Problem. 

^ (Impact  Velocity  = 1 km/s,  Rad  of  Projectile  =0.5  cm) 
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Figure  5.  Compression  Map  at  60/i.s  for  the  Typical  Case. 

(Impact  Velocity  = 1 km/s,  Rad  of  Projectile  = 0.5  cm) 
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Figure  6.  Pressure  Map  at  60  fis  for  the  Typical  Case. 
(Impact  Velocity  = 1 km/s,  Rad  of  Projectile  = 0.5  cm) 
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Figure  9.  Specific  Internal  Energy  Map  at  60 /iS  for  the 
Typical  Case. 

(Impact  Velocity  = 1 km/s,  Rod  of  Projectile  = 0.5  cm) 
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Figure  10.  Map  of  the  Highest  Temperature  of  Cell  at  60  fis 
for  the  Typical  Case. 

(Impact  Velocity  * 1 km/s,  Rad  of  Projectile  * 0.5  cm) 
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,«  Figure  11.  Map  of  the  Highest  Explosive  Temperature  of  Cell  at 

60  (li  for  the  Typical  Case. 

(Impact  Velocity  = 1 km/s,  Rad  of  Projectile  = 0.5  cm) 
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Figure  12,  Plot  of  Temperature  Rise  of  S 
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(1)  The  highest  shock  temperature  in  the  explosive  is  higher  than 
the  highest  shock  temperature  in  the  metal.  ITiis  is  so  because  the 
explosive  is  more  compressible  than  steel. 

(2)  More  plastic  work  is  done  in  the  steel  than  in  the  explosive, 
because  of  the  higher  yield  strength  of  the  steel.  At  later  times,  the 
temperature  in  both  materials  would  have  been  slightly  higher  if  residual 
shock  temperatures  were  taken  into  account. 

(3)  The  highest  plastic  temperature  in  the  explosive  was  asso- 
ciated with  a tracer  which  was  initially  positioned  adjacent  to  the 
metal  interface.  Even  higher  temperatures  would  have  resulted  if 
conduction  from  the  metal  were  taken  into  consideration. 

(4)  Plastic  work  levels  off  in  any  material  either  because  of 
tensile  failure  of  the  material  or  material  melting.  In  the  tensile 
failure  case,  both  the  hydrostatic  and  the  deviatoric  portions  of  the 
stress  in  the  cells  were  zeroed  out.  The  material  near  the  crater  lip 
usually  fails  in  tension.  Melting  of  the  material  is  based  on  its 
specific  internal  energy  acquiring  a value  equal  to  or  greater  than  an 
inputted  quantity  'melt  energy'.  According  to  equation  (9),  yield 
strength  vanishes  when  the  material  melts.  No  additional  plastic  work 
can  therefore  be  done  on  the  material  after  it  melts.  In  our  calculation 
there  is  some  evidence  of  explosive  melting.  The  cell  internal  energies, 
however,  always  remained  lower  than  the  steel  melt  energy;  tensile  failure 
rather  than  melting  was  therefore  the  reason  for  plastic  work  leveling- 
off  in  the  steel. 

Figures  13  through  15  are  plots  of  pressure  versus  time.  In 
Figure  13,  the  pressure  at  the  explosive/steel  interface  is  plotted 
while  in  Figure  14  the  pressure  is  plotted  for  the  interface  between  the 
projectile  and  the  steel  box.  In  about  2 ps,  the  pressure  at  the 
projectile  interface  drops  to  an  insignificant  level  because  of  the  relief 
waves  arriving  from  the  projectile  free  surface  and  the  explosive 
interface.  The  pressure  on  the  explosive/steel  interface,  on  the  other 
hand,  remains  significant  for  much  longer  time.  The  explosive  near  the 
interface  first  fails  at  about  10  ps  but  heals  soon  after.  At  about 
40  ps,  the  crater  lip  is  very  well  formed,  material  near  the  interface 
goes  in  tension  and  fails.  The  pressure  throughout  the  sample  drops  to  an 
insignificant  level  at  this  time,  the  failed  material  therefore  cannot 
heal.  The  plastic  temperature  therefore  does  not  change  after  40  ps. 

The  highest  grid  pressure  vs  time  is  plotted  in  Figure  15.  This  figure 
illustrates  the  dissipation  of  shock  inside  the  target. 

B.  Parametric  Study 

The  test  matrix  shoivn  in  Table  III  was  employed  for  investigating 
the  effects  of  variation  of  impact  velocity,  projectile  radius  and 
casing  thickness.  In  this  study,  the  length  of  the  projectile  was 
always  2.0  cm. 
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i 

\ A _ A 


Calculation 

Casing 

Thickness 

d 

Radius  of 
Projectile 
r 

Impact 

Velocity 

V 

Residual 

Velocity 

V 

RES  ■ 1 + t/Z 

cm 

cm 

km/s 

km/s 

0.25 

I 

0.50 

0.50 

1.00 

0.80 

1.00 

1.00 

0.80 

II 

0.50 

0.50 

2.00 

1.60 

3.00 

2.40 

III 


1.00 


0.50 


1.20 


0,80 


. \ 


1.  Effect  of  Projectile  Size.  Figure  16  shows  the  variation  of 
plastic-work  temperature  with  time  in  the  explosive  and  the  steel  for 
three  projectile  radii,  0.25,  0.50  and  1.00  cm.  Impact  velocities  for 
all  the  three  cases  were  1.0  km/s.  It  is  clear  from  this  plot  that  the 
dependence  of  the  highest  plastic  work  temperature,  as  well  as  its  rise 
time,  on  the  projectile  diameter,  is  weak  and  complicated.  Total 
plastic  work  for  each  of  the  material  packages  is,  on  the  other  hand,  an 
increasing  function  of  the  projectile  radius,  as  seen  from  Figure  17. 

As  the  projectile  radius  increases,  bigger  and  bigger  regions  in  the 
target  are  involved  in  the  deformation  process  and  acquire  high  tem- 
peratures. The  probability  of  getting  ignition  due  to  plastic  work 
therefore  seems  to  be  independent  of  projectile  radius,  as  long  as  the 
fragment  is  large  enough  to  insure  plastic  working  of  the  explosive,  but 
extent  of  the  reaction  may  depend  on  radius.  The  numerical  value  of  the 
highest  temperature  obtained  in  the  explosive  or  the  steel,  however, 
does  not  change  because  the  plastic  work  per  unit  volume  does  not. 

Figure  18  is  a plot  of  the  highest  explosive  shock  temperature 
versus  projectile  radius.  The  calculations  indicate  that  the  shock 
temperature  increases  rapidly  as  the  projectile  radius  increases. 

However,  the  shock  temperatures  must  plateau  when  the  projectile  radius 
is  sufficiently  great  so  that  the  shock  reaching  the  explosive  is  planar. 
In  the  present  calculations,  this  should  have  occurred  at  a radius  of 
about  1 cm,  and  the  failure  of  the  calculations  to  reflect  this  must  be 
due  to  the  finite  zone  size. 

2.  Effect  of  Impact  Velocity.  Figure  19  is  a plot  of  plastic  work 
temperature  in  the  explosive  and  the  steel  versus  time  for  impact 
velocities,  1,0,  2.0,  and  3.0  km/sec.  The  projectile  radius  for  the 
three  cases  was  0.5  cm.  It  is  clear  from  Figure  19  that  the  highest 
plastic-work  temperature  attained  in  either  material  is  not  strongly 
dependent  on  the  projectile  impact  velocity  in  the  range  of  values 
employed  here.  However,  the  plastic  work  temperature  for  the  impact 
velocity  of  3 km/s  leveled  at  a much  lower  value  because  of  an  early 
tensile  failure  of  the  deformational  regime.  Since  the  plastic  work  has 
also  been  seen  not  to  depend  on  the  projectile  radius,  it  can  be  con- 
cluded that  the  highest  plastic  work  temperature  is  independent  of 
projectile  kinetic  energy,  as  long  as  this  energy  is  above  a certain 
minimum  value.  Temperature  ceases  to  rise  after  t = 40  ps  because  of 
material  failure  in  tension.  Figure  20  is  a plot  of  highest  explosive 
shock  temperature  versus  impact  velocity.  The  explosive  shock  tem- 
perature increases  as  the  impact  velocity  increases,  For  a given 
projectile  kinetic  energy,  the  variation  in  shock  temperature  is  more 
dramatic  with  projectile  velocity  than  with  projectile  radius, 

3.  Effect  of  Casing  Thickness.  One  calculation  was  done  using  a 
casing  thickness  of  t = 1.0  cm,  and  impact  velocity,  v = 1.2  km/s.  The 
higher  impact  velocity  was  chosen  so  that  the  projectile  residual 

velocity^ be  the  same  for  the  thicker  casing  as  for  the  typical  case 

11.  R.  F.  Reaht  and  T.  W.  Ipaon,  Trana  ASME,  8R,  364,  1963. 
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Figure  17.  Tot-cil  Plastic  Work  in  Explosive 
vs.  Projectile  Radius  at  t = 40 /xs 
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Figure  19,  Metal  and  Explosive  Plastic  Work  Temperatures 
vs.  Time  for  V = 1,  2,  3 km/s 

(Projectile  Radius  =0.5  cm) 


20.  Explosive  Shock  Temperature 
vs.  Projectile  Velocity 
( Rad.  of  Projectile  = 0.5  cm) 
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with  d = 0.5  cm  and  v = 1.0  km/s.  The  highest  shock  temperature  and  the 
highest  plastic  work  temperature  are  plotted  as  a function  of  the  casing 
thickness  in  Figure  21,  for  both  the  explosive  and  the  casing.  Explo- 
sive shock  temperature  is  seen  to  drop  considerably  as  the  casing 
thickness  increases.  Explosive  plastic-work  temperature,  on  the  other 
hand,  shows  an  increase  with  the  increasing  thickness  which  can  be 
attributed  to  the  delayed  tensile-failure  of  the  explosive  when  thicker 
casing  is  employed.  The  shock  temperature  in  steel  is  slightly  higher 
for  the  thicker  casing,  because  of  the  higher  impact  velocity  while  its 
plastic  work  temperature  is  considerably  higher,  again  because  of  its 
delayed  tensile-failure. 

C.  Effect  of  Reduction  of  Cell  Size 


It  is  generally  believed  that  the  smaller  the  cell  size,  the  greater  is 
the  accuracy  of  the  computational  results.  In  order  to  see  the  effect 
of  reduction  of  cell  size  on  the  numerical  value  of  the  results,  the 
cell  size  employed  in  the  typical  calculation  was  halved  in  the  impor- 
tant deformation  regime.  It  should  be  mentioned  that  the  zone-size 
reduction  is  usually  limited  by  the  limited  computer  core  and  the  rising 
cost  of  computation  (core  size  is  roughly  inversely  proportional  to  the 
square  of  the  cell  size).  Care  must  be  taken  not  to  have  cells  in  the 
grid  whose  aspect  ratio  is  greater  than  two.  If  variable  zoning  is 
used,  the  zone  size  should  change  by  no  more  than  10%  from  one  cell  to 
the  next.  Also,  one  should  have  no  fewer  than  three  cells,  in  order  to 
describe  a particular  material  dimension  in  the  grid.  Figure  22  is  a 
plot  of  the  explosive  and  the  steel  plastic-work  temperatures  versus 
zone  size.  Temperatures  in  both  materials  are  seen  to  increase  as  the 
cell  size  decreases.  It  is  felt  that  there  may  be  a region  near  the 
axis  which  obtains  very  high  strain.  The  code's  ability  to  resolve  such 
a region  is  limited  due  to  the  cell  size  in  the  radial  direction,  which 
has  not  been  altered.  The  cell  size  in  this  calculation  was  halved  by 
shortening  its  axial  dimension. 

If  a linear  extrapolation  is  assumed  to  be  valid,  one  can  easily 
calculate  a plastic  temperature  rise  of  about  100°C  in  the  explosive  and 
about  400°C  in  the  steel,  for  very  small  cell  size.  Since  under  the 
ambient  conditions  the  explosive  melts  at  about  80°C,  a temperature 
higher  than  80°C  will  be  impossible  for  the  explosive  with  the  code 
unless  viscous  effects  are  included.  The  linear  extrapolation  for  the 
explosive  is  therefore  invalid,  however  an  extrapolated  temperature  rise 
of  400°C  for  steel  is  quite  interesting.  It  may  mean  that  a temperature 
rise  of  400'’C,  for  the  solid  explosive  in  the  neighborhood  of  the  hot 
steel,  is  also  possible.  Local  high  temperature  in  steel  will,  however, 
be  shortlived  because  of  its  high  thermal  conductivity. 

D.  Effect  of  High  Tensile  Strength 


The  bulk  strain  for  material  failure  was  decreased  from  -.001  to 
-.1  for  the  explosive  and  from  -.02  to  -.1  for  the  steel,  in  one  of  the 
calculations.  The  change,  which  was  tantamount  to  increasing  the 
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tensile  strength  of  the  material,  did  not  noticeably  alter  the  plastic- 
work  temperature  of  the  explosive.  This  indicates  that  the  shear 
stresses  were  negligible  at  the  time  material  failed  in  tension.  The 
plastic  work  temperature  would,  of  course,  increase  if  the  material 
yield  strength  were  changed,  which  in  turn  indicates  the  significance  of 
work  hardening  terms  of  equation  (9) . 

E.  Effect  of  Presence  of  Cracks  or  Notches  in  the  Explosive  Material 

It  is  well  kno\m  that  crack  tips  act  as  stress-concentrators.  The 
shear  stress  and  the  plastic  work  is  therefore  quite  large  at  the  crack 
tips.  If  the  work  hardening  and  the  effect  of  pressure  on  melting  point 
are  disregarded,  it  follows  that  the  highest  plastic  work  temperature 
attainable  at  the  stationary  crack  tip  in  explosive  is  80°C,  which  is 
the  melting  point  of  the  explosive.  On  the  other  hand,  if  there  is  any 
sliding  between  the  cracked  surfaces  of  the  explosive  material,  it  may 
be  possible  for  the  temperature  to  rise  at  the  cracked  explosive  surface 
by  frictional  work  since  melting  and  heat  conduction  in  the  explosive 
arc  slow. 


IV.  CONCLUSIONS 

Based  on  the  results  obtained  from  the  numerical  simulation  of  the 
normal  impact  of  the  confined  explosives,  one  may  conclude  that: 

(1)  Plastic  work  in  explosive  is  limited  by  melting  and  tensile- 
failure. 

(2)  Plastic  work  is  dependent  on  the  value  chosen  for  the  yield 
strength  of  the  material.  Work-hardening  may  play  an  important  role  in 
determining  the  plastic-work  temperature. 

(3)  If  material  failure  in  tension  is  disregarded,  the  highest 
plastic  work  temperature  in  explosive  or  steel  is  not  too  dependent  on 
the  projectile  velocity  or  the  projectile  radius  values  employed  in  the 
numerical  calculations  here.  For  impact  velocity  of  3 km/s,  however, 
the  material  failed  in  tension  at  a very  early  time. 

(4)  The  highest  plastic  work  temperature  in  the  explosive  leveled 
off  at  about  40  us  regardless  of  the  projectile  velocity  and  the  pro- 
jectile radius.  At  this  time  the  crater  lip  was  well  formed  and  material 
in. the  vicinity  of  the  lip  failed  in  tension. 

(5)  The  highest  shock  temperatures  in  both  steel  and  explosive  are 
strongly  dependent  on  the  projectile  radius  and  velocity. 

(6)  Thicker  confinement  always  lowers  the  explosive  shock  tem- 
perature but  does  not  always  alter  the  plastic  deformation  significantly. 
The  plastic  work  thus  may  become  relatively  more  important  for  the 
thicker  casing. 
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(7)  Explosive  hot  spots  in  the  neighborhood  of  rapidly  deforming 
steel  casing  may  have  a temperature  of  about  400°C.  This  deduction  has 
been  made  by  linearly  extrapolating  the  plastic  temperature  versus  cell 

2 

size  plot  to  the  cell  size  of  0.0  cm  . 

(8)  In  addition  to  plastic  deformation,  heat  conduction  and 
friction  between  the  grits  or  cracked  surfaces  of  the  explosive  may  also 
be  important  heating  mechanisms  in  the  explosives  leading  to  violent 
reactions. 


V.  RECOMMENDATIONS  FOR  FUTURE  WORK 

A parametric  study  of  the  dependence  of  the  explosive  plastic-work 
on  the  projectile  radius  and  the  projectile  velocity  was  attempted  in 
this  effort.  This  study  ignored  the  viscous  contribution  to  the  stress- 
strain  relation.  Addition  of  this  term  might  affect  the  results, 
particularly  the  dependence  of  plastic  work  on  velocity.  In  addition, 
because  of  the  strong  sensitivity  of  the  plastic-work  temperature  on  the 
value  of  the  yield  strength,  it  is  important  to  include  the  work-hardening 
and  the  thermal  softening  terms  in  the  numerical  computational  scheme. 
Other  important  items  that  should  be  considered,  in  any  calculations 
related  to  the  explosive  deformation,  are  the  dependence  of  the  melting 
point  and  the  GrUneisen  constant  on  pressure.  Inclusions  of  failure 
modes  other  than  the  tensile  mode  should  also  be  investigated  and,  if 
needed,  incorporated  in  the  numerical  schemes.  Frictional  work  may  be 
an  additional  important  heating  mechanism  in  the  explosives.  The  slide 
line  concept  employed  in  the  plugging  version  of  HELP  code  may  be  useful 
in  studying  the  relative  material  motion  leading  to  viscous  heating  of  a 
small  layer  of  the  explosive  which  has  just  melted  inside  a narrow  crack 
in  the  explosive. 
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APPENDIX  A.  EQUATION  OF  STATE  OF  COMP  B3 

In  order  to  study  the  shock  wave  propagation  in  materials,  a 
relationship  between  the  equilibrium  thermodynamic  state  parameters,  P, 

V and  E is  required.  Such  a relationship  is  called  "Equation  of  State". 
Many  different  forms  of  these  relationships  are  available  in  literature; 
Gray,  Tillotson  or  Mie-Griineisen  type  equations  of  state  are  more  suit- 
able for  describing  the  thermodynamical  state  of  a solid  explosive,  such 
as  Comp  B3. 

In  this  report  an  equation  of  the  Mie-Griineisen  type  has  been  used. 
This  equation  was  calibrated  by  employing  an  experimental  Hugoniot.  In 
addition,  the  Griineisen  parameter  was  assumed  to  be  a constant.  A 
Hugoniot  typically  defines  all  the  pressure-volume  states  obtainable 
through  a shock  transition.  Usually  the  Hugoniot  of  material  is  avail- 
able in  the  form  of  experimental  data  points  consisting  of  shock  velocity 
and  particle  velocity.  Usage  of  Rankine-Hugoniot  jump  conditions 
reduces  the  above  data  to  equivalent  data  consisting  of  pressure  and 
volume  (or  compression).  Following  sections  illustrate  the  derivation 
of  a Mie-Griineisen  equation  of  state  for  Comp  B3  starting  from  the  shock 
wave  data. 

A.  Shock  Wave  Data 


N,  L.  Coleburn  and  T.  P.  Liddiard’^  first  reported  the  shock  wave 
data  for  Comp  B3  which  is  reproduced  in  Table  A-I. 

This  experimental  data  was  subjected  to  a least  square  curve  fit 
yielding  the  following  relationship  between  the  shock  velocity  and  the 
particle  velocity. 


Ug  = 2.66  + 1.92  U^^  (km/s)  (A-1) 

Equation  (A-1)  expresses  the  material  shock  velocity  as  a linear  function 
of  the  particle  velocity.  The  constant  term  in  equation  (A-1)  has  the 
significance  of  sound  speed,  c^.  Equation  (A-1)  also  expresses  the  fact 

that  the  sound  speed  is  the  lowest  value  a shock  can  acquire.  Ideally  a 
material  subjected  to  a plane  shock  undergoes  no  lateral  deformation. 

Each  macroscopic  volume  element  is  compressed  exactly  as  its  lateral 
neighbors  and  is  therefore  deformed  only  in  the  direction  of  the  shock 
propagation.  In  reality  the  material  develops  a shear  and  if  the  yield 
stress  is  exceeded,  plastic  flow  takes  place  in  the  material  on  a micro- 
scopic scale.  As  a result,  the  stress  of  the  material  cannot  just 
consist  of  isotropic  hydrostatic  pressure,  but  will  have  a shear  com- 
ponent as  well.  In  the  experiments,  however,  only  the  normal  component 
of  the  stress  is  usually  measured,  the  shear  has  to  inferred  from  the 
material  constitutive  relationship. 


* N.  L.  Colehum  and  T.  P.  Liddiavdy  J.  Chem.  Phys.y  4£,  1929,  2996. 
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le  A-tI.  Shock  Wave  Da 

ta  for  Comp  B3  (p^  = 

km/s 

km/s 

3.387 

0.380 

3.437 

0.421 

3.439 

0.469 

3.510 

0.466 

3.628 

0.462 

3.713 

0.460 

3.971 

0.741 

4.243 

0.819 

4.314 

0.886 

4.346 

0.869 

4.377 

0.884 

4.441 

0.936 

4.469 

0.881 

The  basic  equations  used  in  the  study  or  shock  are  obtained  by 
straightforward  application  of  the  conservation  of  mass,  momenta  and 
energy  across  a shock  front.  Assumption  is  made  here  that  the  initial 
and  the  final  states  are  equilibrium  states.  In  Figure  A-1,  a shock 
moving  into  the  material  at  rest,  with  a velocity,  U^,  produces  a 

particle  velocity,  U^,  and  raises  the  pressure  from  to  Pj^,  the 

density  from  to  Pj^,  and  the  specific  internal  energy  from  to  Ej^j. 

The  Rankine-Hugoniot  relationships  expressing  the  conservation  laws 

(A-2) 


are 


P U 
0 s 


u ) 

p 


P„  - P = p u u 
H 0 *^0  s p 


(A-3) 


- ^0  = i ^ ^0^  ("h 


V ) 


(A-4) 


Equations  (A-2)  and  (A-3)  can  be  rearranged  to  yield  expressions  for 
pressure  and  compression, 


p U U 
o s p 


(A-5) 


""  p 


- 1 = 


U - U 
s p 


(A-6) 


where  the  initial  pressure  P^  has  been  neglected.  Equation  (A-1)  can 
now  be  used  to  eliminate  Up  from  equations  (A-5)  and  (A-6),  yielding 


Pll  = 


P U (U 
o s ^ s 


2.66) 


1.92 


(A-7) 


and 


U - 2.66 
s 

^ 0.92  U + 2.66 


(A-8) 


For  a given  shock  speed,  equations  (A-7)  and  (A-8)  yield  a data  point. 

A number  of  such  data  points  are  obtained  for  shock  velocities  in  the 
range  ^ 3c^.  These  values  can  be  fitted  to  a polynomial  of  the 

form 


Ph  = Fp^  + 


(A-9) 
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Figure  Al.  A Plane,  One -Dimensional  Shock  Wave 
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Equation  (A-9)  gives  the  value  of  pressure  for  any  Hugoniot  shock 
compression  state,  Once  the  pressure,  P^j,  is  determined  from 

equation  (A-9) , the  Hugoniot  energy  can  simply  be  obtained  from 
equation  (A-4) , 


B.  Mie-Grlineisen  Equation  of  State 

A Hugoniot  state  characterized  by  Ppj  (or  y^j)  and  E^j  can  now  be  used 

as  a ready  reference  for  determining  all  the  possible  material  states. 
These  states  may  not  necessarily  be  reached  by  a plane  shock  compression. 
Such  states  may  arise  in  a real  material,  where  the  finite  dimensions  of 
the  sample  will  give  rise  to  rarefactions  or  in  the  case  where  the 
impacting  projectile  and  the  resulting  shock  are  non-plan.ir.  In  these 
cases,  lateral  deformations  are  not  negligible.  A convenient  equation 
of  state  is  of  the  Mie-Grtlneisen  type: 

P - Ph  = (E  - (A- 10) 


where  F(V)  is  a Grlineisen  parameter,  P,  V and  E are  any  equilibrium 
pressure,  volume  and  energy.  An  expression  for  the  Grlineisen  parameter 
is  obtained  by  differentiating  equation  (A-10), 


r(V)  = V (f)  = ^ (§) 

V V V 


(A-11) 


Usually  the  Grlineisen  parameter  is  a weak  function  of  volume,  its  volume 
dependence  will,  therefore,  be  ignored.  On  substituting  equation  (A-4) 
into  (A-10),  one  obtains 


P = Ph  (1  - ^)  + rpE  (A-12) 

where  p = 1/V  is  the  density  of  the  material  corresponding  to  pressure  P 
and  energy  E.  On  substituting  the  expression  for  Ph  from  the  polynomial 

fit  of  equation  (A-9)  into  (A-12),  one  obtains 

P = (Ay  + By^  + Cy^)  + PpE  (A-13) 

where  A,  B,  C are  constants  given  by: 

A = 13.5  X 10^° 

B = 9,5  X 10^°  (dyno/cm^)  (A-14) 

C = 100.6  X 10^® 
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Equation  (A-13)  is  a general  equation  of  state  used  to  describe  the 
response  of  Comp  B3  under  dynamic  loading.  Figure  A-2  is  a schematic 
showing  the  relationship  between  the  equilibrium  thermodynamic  parameters, 
P,  V ajid  E.  In  Figure  A-2,  a shock  hydrostat  is  shown  along  with  a few 
pressure-volume  curves  with  constant  energies. 

It  should  be  pointed  out  that  equation  (A-13)  is  valid  strictly  for 
the  compressed  states,  i.e.  for  y > 0.  The  y < 0 values  correspond  to 
expanded  states  of  material,  for  which  the  hydrostatic  pressure  must 
become  tensile.  In  order  to  ensure  this  behavior  from  an  equation  of 
state,  it  is  customary  to  change  the  sign  of  the  square  term  in  equa- 
tion (A-13).  The  resulting  curve  is  shown  in  Figure  A-3.  Comp  B3  fails 
for  small  negative  pressure.  Appendix  B indicates  that  Comp  B3  fails  for 
y £ -.001. 

It  is  generally  assumed  that  a Mie-Grlineisen  type  equation  of  state 
is  adequate  for  metals  in  the  low  pressure  regime,  P < 10  kbar.  Typical 
peak  pressures  obtained  in  impacts  involving  confined  explosives  are 
usually  less  than  50  kbar.  The  suitability  of  Mie-GrUnoisen  equation  to 
describe  the  explosive  behavior  is,  however,  still  an  open  question. 


APPENDIX  B.  MATERIAL  CONSTANTS  FOR  COMP  B3 


This  section  lists  the  values  of  miscellaneous  parameters  that  are 
usually  needed  as  input  for  a computational  work  involving  Comp  B3. 
These  values  have  been  collected  from  several  sources.  In  the  present 
computational  work,  the  parameters  of  the  explosive,  PBX  9404  were 
employed  if  the  corresponding  parameters  for  Comp  B3  were  not  readily 
available. 

MISC  CONSTANTS 


3 

Density,  = 1.68  g/cm 
Coefficient  of  friction,  ii  = .35 

2 

P t 2 

Critical  shock  initiation  energy,  — rj—  = 29  cal/cm 

P 


SOUND  VELOCITIES 


Longitudinal  shear  velocity 


^ r 
3 ^ 


3,15  X 10^  cm/s 


Transverse  shear  velocity, 


Bulk  sound  speed, 


C = J-  = 1.73  X 10^  cm/s 
sip 

= 2.44  X 10^  cm/s 


ELASTIC  CONSTANTS 

2 10  2 

Bulk  Modulus,  ^ = Pq  10.037  x 10  dyne/cm 

2 10  2 

Shear  Modulus,  G = p^  c^  = 5.029  x 10  dyne/cm 

Poisson's  ratio,  v = . = .29 

2(G  + 3tc) 

Hugoniot  elastic  limit,  = 1.5  kb  (PBX  9404) 

Tensile  yield  stress,  = (|~^)  = 0.88  kb 

Shear  yield  stress,  Y„  = — ^ Y„  = 0.51  kb 

^ 3 ' 


Ultimate  tensile  strength,  S = .1042  kb 

S 

Volumetric  strain  for  failure,  " " ic  “ 

Maximum  distension  for  failure,  AMDM  = + 1 = .999 

THERMAL  PROPERTIES 

Normal  melting  temperature, 

Energy  required  to  melt,  E^  = J '^MELT  ~ ^ ®rg/g 

30°C 
50°C 
70°C 

83°-100°C 

Thermal  conductivity,  k = 6.27  x 10  ^ cal/cm-sec-°C 

Coefficient  of  thermal  expansion,  a = 54,6  - 97.5  x 10  cm/cm-0°C 

9P 

GrUneisen  constant,  r = C'^)  = 0.947 


Specific  heat,  = 


.299  at 
.307  at 
.325  at 
\ .333  at 


APPENDIX  C.  WALSH- CHRISTIAN  TECHNIQUE 


Temperature  rise  due  to  shock  compression  of  metals  along  a Hugoniot 
and  subsequent  temperature  drop  of  the  material  along  an  adiabat  can  be 
easily  calculated  by  simple  application  of  thermodynamics.  The  thermo- 
dynamical identity. 


T dS  = C^  dT  + T dV 

V 


(C-1) 


The  Rankine-Hugoniot  jump  condition. 


^H  - ^o  = i (>^H  ^ '’o^  (V 

and  the  second  law  of  thermodynamics. 


tC-2) 


T dS  = dE  + P dV 


(C-3) 


form  a starting  point  for  such  a calculation.  The  derivation  of  an 

appropriate  expression  was  given  by  Walsh  and  Christian^^  several  years 
ago  and  is  presented  here  for  the  sake  of  completeness. 


From  equations  (C-1),  (C-2)  and  (C-3),  it  is  easily  seen  that  the 
energy,  along  a Hugoniot,  is  not  an  independent  variable  and  can 
therefore  be  eliminated.  This  is  accomplished  by  first  integrating 
equation  (C-3)  and  then  substituting  equation  (C-2)  for  the  energy 
terms.  The  result  is 


^ 'V  ‘'^Ihug  = I <-h  * V - V * { 
^0  '^0 


(C-4) 


Integration  of  equation  (C-1)  yields 


/"(T  dS]  = C„  (T„  - T„)  * T dV 

S V V 

0 0 


(C-5) 


Differentiating  equations  (C-4)  and  (C-5)  and  equating  the  right  hand 
sides  of  equations  (C-4)  and  (C-5)  eliminates  the  entropy  terms  as  well. 
The  result  is 


where 


d T„  8 P,, 

^ (dir)  ^ \ ^ 

H 


(C-6) 
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(C-7) 


(V)  = I P 


- 

2 dV  '■'^0 


V) 


Further  simplification  of  equation  (C-6)  results  by  assuming  that  C and 
9P 

(■^)  are  constants.  The  resulting  equation  is  a nonhomogeneous,  first 

V 

order  differential  equation 


d T„  f CV„) 

“ - b T„  = " 


d V„  ‘H  C 
H \ 


where 


b = 


(9P/9T) 


The  solution  of  equation  (C-8)  is  easily  obtained  to  be 


(C-8) 


(C-9) 


- V„)  -bV,,  V„  bV 

T„  = T.  e ° 


H 


V 


HUG 


CC-10) 


The  second  term  on  the  right  hand  side  of  equation  (C-10)  can  be  expressed 
in  terms  of  Hugoniot  parameters;  A,  B,  C [Equation  (1)].  The  result  is 


T 


H 


+ 


1 _£ 

2 C 


bVo 

.‘'h  (A  + B)u^  + 2(B  + C)p^  + 3C|j'’  Uu  , 

j g (jj] 

0 (1  + V) 


(C-11) 


The  integral  in  equation  (C-11)  can  be  evaluated  analytically  but  is 
obtained  here  numerically  in  order  to  save  time.  It  should  be  mentioned 
that  in  this  work,  only  the  peak  temperature  obtained  by  explosive 
during  the  shock  is  of  concern.  However,  soon  after  the  passage  of  the 
shock  wave,  material  begins  to  release  itself  along  an  adiabat.  The 
temperature  of  the  material  as  a consequence  of  this  adiabatic  expansion, 
drops.  The  temperature  corresponding  to  any  volume  attained  by  the 
material  during  this  expansion  is  easily  calculated  from  equation  (C-1), 
where  dS  is  set  equal  to  zero.  The  resulting  equation  is 


C dT  + 

V 


T dV  = 0 


(C-12) 
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The  solution  to  this  is  given  by 

''ad  ^ 

Normally,  is  the  inverse  of  density  of  the  material  under  full 

compaction.  Equations  (C-11)  and  (C-13)  give  a complete  temperature 
history  of  the  material  undergoing  a shock  compression  and  subsequent 
adiabatic  expansion.  Because  of  the  irreversible  nature  of  the  shock, 
both  the  energy  and  temperature  of  the  system  are  higher  after  the 
passage  of  the  shock,  the  residual  temperature  being  given  by: 


T = T 
RES  AD 


(C-14) 
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LIST  OF  SYMBOLS 


Symbol 


Definition 


Units* 


AMDM 

b 


d 

E 

E. 


J 

P 

P. 


^ij 


compression  for  tensile  failure  of  material 

r/v 

sound  speed 

specific  heat  at  constant  volume 
casing  thickness 
specific  internal  energy 
Hugoniot  specific  internal  energy 
melt  energy  per  unit  mass 
mechanical  equivalent  of  heat 
pressure 

Hugoniot  pressure 
projectile  radius 
deviatoric  stress 
spall  strength 
time 

temperature 
initial  temperature 


-3 


cm 
km/s 

cal/g-°C 

cm 

erg/g 

erg/g 

erg/g 

erg/cal 

dyne/cm^ 

dyne/cm^ 

cm 

dyne/cm^ 

dyne/cm^ 

S 


* The  CGS  system  of  units  is  employed  here  in  order  to  he  consistent 
with  the  input/output  of  the  HELP  code. 
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Symbol 

Definition 

Units 

u. 

] 

velocity 

cm/s 

U 

P 

particle  velocity  behind  the  shock 

km/s 

U 

s 

shock  velocity 

km/s 

V 

impact  velocity  of  the  projectile 

km/s 

Res 

residual  projectile  velocity 

km/s 

V 

specific  volume 

3 

cm 

initial  specific  volume 

3 

cm 

W 

specific  plastic  work 

erg/g 

Y 

yield  strength 

dyne/cm 

r 

Grllneisen  constant 

-- 

e.  . 

strain 

-- 

n 

p/Pq 

-- 

K 

bulk  modulus 

dyne/cm' 

0.  . 
IJ 

stress 

dyne/cm' 
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